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ABSTRACT 

In the absence of an active dynamo, purely poloidal magnetic field configurations are 
unstable to large-scale dynamical perturbations, and decay via reconnection on an 
Alfvenic timescale. Nevertheless, a number of classes of dynamo-free stars do exhibit 
significant, long-lived, surface magnetic fields. Numerical simulations suggest that the 
large-scale poloidal field in these systems is stabilized by a toroidal component of the 
field in the stellar interior. Using the principle of conservation of total helicity, we 
develop a variational principle for computing the structure of the magnetic field in- 
side a conducting sphere surrounded by an insulating vacuum. We show that, for a 
fixed total helicity, the minimum energy state corresponds to a force-free configura- 
tion. We find a simple class of axisymmetric solutions, parametrized by angular and 
radial quantum numbers. However, these solutions have a discontinuity in the toroidal 
magnetic field at the stellar surface which will exert a toroidal stress on the surface 
of the star. We then describe two other classes of solutions, the standard spheromak 
solutions and ones with fixed surface magnetic fields, the latter being relevant for neu- 
tron stars with rigid crusts. We discuss the implications of our results for the structure 
of neutron star magnetic fields, the decay of fields, and the origin of variability and 
outbursts in magnetars. 
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1 INTRODUCTION 

The importance of magnetic helicity has long been ap- 
preciated in the c ontext of force-free plasma dynamics 
(|Woltierlll9"53 . ll962l 'l. This is due primarily to the fact that 
helicity has been experimenta lly shown to be conserved 
during magnetic reconn e ction |Tavloi Il974l: |ji et al.l 1 19951 : 
iHeidbrink fc Dandbood : IHsu fc Bellanll2002l) . even when a 
substantial fraction of the magnetic energy has been dissi- 
pated. 

In the context of stars, numerical magnetohydrody- 
namics simulations of initially random fields in a con- 
ducting sphere reveal that magnetic helicity plays a cru- 
cial role in d etermining the final magnetic configuration 
jBraithwaite fc Spruit 2004). In particular, after a period 
of violent reconnection, which occurs on the Alfven cross- 
ing timescale, a stable, mostly dipolar field develops, the 
strength of which depends solely upon the initial helicity. 
In the absence of initial helicity, the final magnetic field 
vanishes. Relevant systems include magnetic A stars, white 
dwarfs and neutron stars. In the present paper we are inter- 
ested primarily in neutron stars. 

The simplicity of the final solution found by 
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Braithwaite fc Spruij l|2004 ) and others (e.g.. lYoshida et al.l 
20061 ) suggests that a simple analytical model can be devel- 
oped. We present here a formalism that might be relevant 
for understanding the nature of the final equilibrium. Our 
approach is based on a variational approach that encodes 
th e principle o f conservation of total helicity, as proposed 
bv lTavloij|l974l . We do not address the dynamical path that 
a system will follow starting from its initial state, but sim- 
ply attempt to construct the final state of the system. Thus, 
we present a set of magnetic field solutions as a function of 
the magnetic helicity. We obtain these solutions via a vari- 
ational energy principle in which the integrated helicity is 
held fixed. 

In 31] we state the variational principle and derive the 
governing field equations. ^ which is the heart of this paper, 
then solves the equations under different conditions, assum- 
ing axisymmetry, and discusses the long term evolution of 
the resulting magnetic field. Explicit forms for some illustra- 
tive solutions are given in ^3.61 ^discusses the implications 
of the results and considers possible applications to neutron 
stars. A number of technical details are given in the Ap- 
pendices, including a comparison of our approach, which is 
based on conservation of total helicity, and the approach of 
citepFiel:86, who considers a conserved microscopic helicity 
associated with individual magnetic flux tubes. 
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2 HELICITY, MINIMUM ENERGY, AND THE 
FORCE-FREE CONDITION 

The magnetic helicity density is defined by ft = A • B, where 
A is the vector potential and B = VxA is the magnetic 
field. Generally, h is gauge dependent. Nevertheless, the vol- 
ume integrated helicity can be made gauge independent by 
a suitable choice of integration. That is, 

i/= / -^A-Bd^a;, (1) 

Jv Stt 

is gauge invariant if B ■ dS = everywhere on the boundary 
of V. While during reconnection the local helicity density is 
not conserve d, to a good approximation the total helicity is 
l|Tavloij |l974:). This experimental fact suggests that the in- 
tegrated quantity H plays an important role in determining 
the final magnetic field configuration, and is the foundation 
of our analysis. In contrast to the helicity, the magnetic en- 
ergy. 



U : 



B B d^'a; , 



(2) 



is not conserved during reconnection. Indeed, reconnection 
can convert a substantial fraction of the magnetic energy 
into thermal energy of the plasma. 

In the absence of a dynamo, the equilibrium state of the 
field will be associated with its minimum energy state. That 
is, reconnection will proceed until the maximum amount of 
magnetic energy is converted into heat. If there are no con- 
straints on the system, the entirety of the magnetic energy 
will be converted. However, even under randomly chosen ini- 
tial conditions, a non-zero helicity will initially be present, 
and subsequently must be conserved. Thus, there is a limit 
to how much magnetic energy can be converted into heat, 
and the final configuration must consist of a non-vanishing 
magnetic field. This situation naturally lends itself to a vari- 
ational approach, which we explore in the following subsec- 
tions. 

In our analysis we seek to minimize the total magnetic 
energy ([2)| subject to a fixed value of the total helicity 
However, before doing this we might ask: Is it correct to min- 
imize only the magnetic energy and to exclude contributions 
to the energy from the underlying fluid? The answer is yes. 
In the context of a constrained total helicity, the equations 
governing the fluid equilibrium, up to a non-radial surface 
stress (which we will assume is balanced by a rigid crust), 
are completely separable from those governing the magnetic 
fleld structure. (The reader is referred to Appendix |X] for 
a discussion of this question and a comparison of our total 
helicity conservation approach and that of Field (1986) who 
considered helicity conservation on individual fleld loops). 
That is, subject to the constraint of flxed helicity, we might 
have considered minimizing 



Uu 



1' 




d\-+ / 


Iv 




Jv 



Bd^x 



Stt (3) 
= Us{p,V) + U{B,V), 

where e and $ are the speciflc internal energy and self- 
gravitational potential of the fluid, respectively, and we as- 
sume for simplicity that the star is barotropic (although 
this is not necessary). Since the helicity constraint depends 
upon only the magnetic fleld, the magnetic and hydrody- 
namic contributions to the energy are distinct, and thus, 



for a fixed V , may be minimized independently. As a conse- 
quence, the equilibrium magnetic field configuration will not 
explicitly depend upon the hydrostatic structure of the star. 
Furthermore, such a state will correspond to a global min- 
imum of the combined energy, Utot- Perturbations around 
such a state will necessarily increase Uj and U , otherwise 
these would not be independent minima, and thus will in- 
crease Utot - Here and henceforth, when we refer to "helicity" 
we mean the total helicity defined in equation ([T|, as dis- 
tinct from the local helicity density ft = A ■ B or the helicity 
associated with a single magnetic loop (Appendix IX)) . 

A somewhat subtle point to note is that the fluid and 
magnetic fleld structure are implicitly coupled by the defl- 
nition of V , the volume of the conducting region (the star). 
That V is flnite is critical to the existence of a stable mag- 
netic fleld conflguration. Indeed, if V is inflnite there can 
exist no equilibrium magnetic fleld conflguration with non- 
vanishing B. This is a result of the fact that U is propor- 
tional to , and hence the field will naturally expand to fill 
all available space and therefore have vanishing strength. In 
the case of a star, equilibrium is a consequence of the star's 
self-gravity. This enters into the determination of the mag- 
netic field by constraining the size of the conducting region, 
V, which anchors the magnetic field. Generally, Utot must be 
minimized with respect to V as well as p and B. However, 
since both U and Uf are functions of V , the stellar radius, 
-R, implicitly couples the structure of the magnetic field and 
the fluid. We will discuss how this may be done explicitly in 

Woltierl |l95^) demonstrated that, in the context of uni- 



form media, extremizing the magnetic energy subject to a 
flxed total helicity results in the force-free equation, 

VxB = qB, (4) 

where a is a Lagrange multiplier (Appendix|Bj. Such a fleld 
configuration is known as a Woltjer state. The situation is 
more complicated if there are constraints upon the currents 
that can be driven in the conducting material. The simplest 
case, and the one of particular interest here, is when the 
system under consideration can be separated into two dis- 
tinct regions: an inner region of infinite conductivity (type-I 
region) which corresponds to the interior of the star, and 
an outer region of vanishing conductivity (type-II region) 
which corresponds to the vacuum exterior of the star. 

Given two such regions, the standard variational prin- 
ciple for uniform media may be extended by the inclusion of 
Lagrange multipliers a and A(a;). Thus, we seek to extremize 
the quantity 

B-B-aA-B + — /A-J)d^a;, (5) 



S ■- 



1 

V 8^ 



where the current J is related to the magnetic field in the 
normal way, J = cVxB/47r, and f{x) vanishes in regions of 
infinite conductivity (type-I region) and is unity otherwise 
(type-II region). Unlike a, which corresponds to fixing H 
and is a scalar constant, the A are vectors and functions of 
position; the difference arises because of the global versus 
local nature of their respective constraints. 

Varying 5* with respect to A, holding SA flxed on some 
outer boundary {not on the boundaries between type-I and 
type-II regions) gives 



VxB = QfB - VxVx/A, 



(6) 
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and varying with respect to A trivially gives 
/VxB = 0. 



(7) 



Within type-I (/ — 0) regions, equations (|6} & (0 reduce 
identically to the Woltjer force-free condition Q. In con- 
trast, within type-II (/ = 1) regions, the additional terms 
are required to make equations ((6)1 and consistent with 
each other. In particular, inserting the latter into the former 



This implies that 



aB = Vx VxA . 



VxA = QfA + aVA 



(8) 



(9) 



for some A, the explicit form of which depends upon the 
gauge choices for A and A. A natural gauge for the former 
is the so-called Force-Free Gauge in which B = aA (see 
Appendix|B]). This happens also to be a Lorentz gauge (i.e., 
since we are concerned with stationary solutions, V-A = 0). 

As discussed in detail in Appendix[Cl the discontinuous 
nature of / at the boundaries between type-I and type-II 
regions implies the standard set of boundary conditions: 



n- (Bi - Bii) =0, 
n X (Bi - Bii) = —K 



(10) 



Generally, for each choice of a and A, a family of solutions 
which are local minima of magnetic energy subject to the 
helicity constraint may be produced. The global minimum 
can then be found by inspection. 

In the case of a uniform medium it is possible to show 
that the magnetic energy and helicity are related by 



U = aH. 



(11) 



In contrast, for non- uniform media no such simple relation 
exists. This is because the force-free gauge is not gener- 
ally a Lorentz gauge in the type-II regions. However, we 
can parametrize the departure from the uniform-medium 
expression in terms of the gauge transformation that relates 
a Lorentz gauge in the type-II regions to that determined 
by the continuity of A at the boundaries. 

In particular, let ip be the field potential of the magnetic 
field in type-II regions (i.e., B — Vi^s) and Al be a gauge 
function which relates A to the Lorentz gauge in type-II 
regions (i.e., A — Al -I- VAl where V-Al ~ 0). Then, as 
shown in Appendix |D] U and H are related by 



U -aH - 



a(i-ii) 



Stt 



(AlVv3-v3VAl) -dS, (12) 



which may be identified as the action. 



3 AXISYMMETRIC SOLUTION 

The discussion in the preceeding section was fairly general. 
However, in the context of stellar magnetic field evolution, 
the case of interest is the magnetic field of a spherical con- 
ducting region surrounded by an insulating vacuum. There- 
fore, henceforth we restrict our attention to the case in which 



if |a;| iJ, 

1 otherwise , 



(13) 



where R is the stellar radius. We explicitly construct the 
solution using vector spherical harmonics (Appendix [E|l: 



Y;„ = erYhn , *i„ EE rVYi 



i 



r X VYir. 



(14) 



We also make the further simplification of assuming that the 
final solution is axisymmetric. 

Assuming axisymmetry simplifies the computation of 
the force-free solutions considerably, since the magnetic field 
can be written as 

B = Vx (yl-^e^) +5-^6^. (15) 

The general procedure consists of solving for A'^ and B'* for a 
force-free magnetic field in the stellar interior and matching 
this to a vacuum field in the exterior. 



3.1 General Solution 

3.1.1 Vacuum Exterior 

In the exterior, the magnetic field is defined by equations 
(fT5l) and 0: 

VxB = VxVx (^A'^e^) + Vx (^B'^e^) = 0. (16) 

We expand A'^e^ and B^e^f, in terms of the vector spherical 
harmonics, 

I I 

where we have dropped the azimuthal quantum number m 
because of axisymmetry (m = 0). We then obtain 



/^lAi = , S; = , 



(18) 



where A; is the 3-dimensional Laplacian associated with the 

(19) 



A; = —OrT Or 5 



This has the general solution 



(20) 



Regularity at infinity requires 6; = 0, and thus the exterior 
field is given by 



B 



EH(' + l)^Y,+/-i|,*,] . (21) 



Note that the exterior B is purely poloidal. 



3.1.2 Stellar Interior 

In the stellar interior it is possible to make use of the force- 
free gauge (B — qA), and rewrite B in terms of B** alone: 



B = a 



-^Vx (S^e^) +B^e^ 



(22) 



The vector field B'^e^ is again expanded in vector spherical 
harmonics, with the result that the force-free equation is 



(23) 



Upon making use of the properties of the vector spherical 
harmonics, this gives 



AiBi+aBi = 0. 



(24) 
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The general solution of equation (|24|l can be written in 
terms of spherical Bessel functions: 



S; = Q [cini{ar) + diji{ar)] 



(25) 



where ji and ni are the spherical Bessel functions of the first 
and second kind, respectively. Regularity at the centre and 
the orthogonality of the n; imply that the ci must vanish. 
Therefore, the general solution for the axisymmetric force- 
free magnetic field in the stellar interior is 



B 



E 



1(1- 



-diji{ar) 



-drrdiji{ar) 



+adiji{ar)^i} . (26) 



In general, the interior field has both poloidal and toroidal 
components. 



3.1.3 Matching Condition 

Equations (|21|l and (|26|) give the general form of the mag- 
netic field in the vacuum exterior and the conducting stellar 
interior. Each is expressed in terms of a set of coefficients, 
ai and di, respectively. The two solutions must be matched 
across the stellar surface, r = R, using the standard bound- 
ary conditions (equation llOf) . This gives the following rela- 
tion between the ai and df. 

0.1 



di = 



(27) 



3.2 Minimum Energy Solutions 

The set of field solutions (|21|l and (I26p with the matching 
condition (|27[) all satisfy the equations obtained from vary- 
ing S with respect to A. However, we have yet to determine 
the minimum energy state. That is, we have replaced the 
degrees of freedom associated with the functional form of B 
with an infinite, but countable, set of degrees of freedom as- 
sociated with the coefficients ai. Furthermore, we have yet to 
choose an a. This can be addressed by considering S{ai; a) 
explicitly, i.e., minimizing the energy with respect to the ai 
and a, subject to the constraint of fixed H. 

Equation (|12p provides a direct way in which to con- 
struct S{ai;a) in terms of (p and Al, and thus we will de- 
termine functional forms for these now. By inspection of 
equation (^1]) . 



E 



ail 



Yin 



(28) 



Determining Al is more difficult. Letting A — A'^etf, is suffi- 
cient to produce the required vacuum field solution, and also 
satisfies the Lorentz gauge. However, A must be continuous 
across the stellar surface. Thus, Al is defined by 



where dzzji{z) 



l{l + l] 



■jiiaR) 



R 



ji{aR) - aji^i{aR) 



(29) 



zji-i{z) — lji{z) was used 



(Abramowitz fc StegunI \l97i ). Setting Al 
j:^Xu{r)Yio{e,^) gives 



VAl=E 



r 



(30) 



which, together with the orthogonality of the vector spheri- 
cal harmonics, implies 



At.! — — 



-(i+i) 



IjiiaR) - aRji_r{aR) 



(31) 



Making use of the above results, and after integrating 
over the stellar surface, we obtain 



5(a..)^£E 



^^1(1 + 1) jj^^ (32) 



R^' ji{aR) 



As expected, variations of S{ai;a) with respect to a gives 
back H{ai;a) (cf. with equation [39]). Variations with respect 
to ai gives a new condition, 

ji-ijaR) 

aai . = , (33) 

ji[aR) 

which implies that either a; = or ji^i{aR) — for all I. 
Since the zeros of the spherical Bessel functions are distinct, 
at most one multipole component can be non-vanishing. 

Thus, we conclude that a minimum energy configura- 
tion consists of a single multipole / with the value of a se- 
lected such that 

ji^i{aR)^0, Oi/O. (34) 

It is worth noting that the above condition on a corresponds 
to the poloidal component of the field being continuous 
across the stellar surface. That is, perhaps unsurprisingly, 
the minimum energy is obtained when there is no kink in 
the poloidal component of the magnetic field at the surface. 
There is, however, generally a kink in the toroidal compo- 
nent of the field. 

The toroidal kink at the surface will exert a toroidal 
stress upon the surface of the star, the direction of which 
depends upon the sense of the magnetic field. In the case 
of a purely fiuid star, the stress will cause motions that will 
result in the untwisting of field lines which cross the surface 
and thus a reduction in the overall helicity of the system. 
One way to eliminate this effect is to focus on solutions with 
zero stress at the surface. We describe a class of such solu- 
tions in i]3.5l Alternatively, if the star has a solid crust, it 
can balance a non-vanishing surface stress by rigidly cou- 
pling different parts of the surface. We discuss the effect of 
a crust in i|3.7l below. For now we continue our discussion 
of the minimum energy solutions given by equation p4p . ig- 
noring the stress at the surface. As it happens, many of the 
qualitative features of these solutions carry through to the 
other situations. 



3.3 Energy and Helicity 

For the values of a; and a found in the previous subsection, 
S identically vanishes (substitute eq.[33]in eq.|32}. As a con- 
sequence, in an energy minimum state, the magnetic helicity 
and energy are related by equation (llip . exactly as for Wolt- 
jer states in uniform media. After integrating equation ([Dip 
for the magnetic energy, we find that 



U = aH 



■ 1(1 + l)a^a^ I 



iTTR^ 



(35) 
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Figure 1. Meridional slices of illustrative magnetic field geometries for (a) the global minimum energy configuration, which corresponds 
to the first energy minimum of a dipolc configuration, (b) & (c) the second and third energy minima of a dipole configuration, and (d) 
the lowest energy minimum of a quadrupole configuration. The poloidal and toroidal field structure are shown by the field lines and 
fiUed-color contours, respectively. In all plots the field line density is proportional to the poloidal flux and is normalized such that the 
average surface fields are equal. Similarly, the toroidal color scheme is uniform among all of the plots. 



For a given fixed helicity H and given I and a, tliis relation 
determines the value of a;. 

We have seen in the previous subsection that we have 
an infinite number of minimum energy solutions. This is be- 
cause we are free to choose any value of the multipole index 
I, and for each I, we may choose any zero of (see ea. l34p 
to evaluate a. Every one of these solutions is a local energy 
minimum with respect to variations of all quantities (includ- 



ing the fluid variables, see §1). Of these many local minima, 
one corresponds to the global energy minimum for the given 
magnetic helicity. Since U = aH, the global minimum will 
have the smallest value of a among all the solutions, i.e., the 
smallest root of ji-i{z) for 1^0. By inspection it is clear 
that the global minimum energy solution has I — 1 and 
a = tt/R, which corresponds to the first zero of jo{z). The 
magnetic field configuration of this solution corresponds to 
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a dipole and is shown in Figure [T^. As mentioned, there are 
countless higher energy solutions, each of which corresponds 
to a local minimum. These involve higher multipoles and/or 
more complex radial structure. Some examples are shown in 
panels b-d of Figure [TJ and are discussed in more detail in 
the following subsection. 

As an aside, we note that these equilibrium configura- 
tions describe a class of solutions, defined by aR and H, as a 
function of R. Since aR and H are fixed, U ocR-^. There- 
fore, if we were free to minimize U with respect to R, we 
would find that R ^ oo and B 0. This does not happen 
in practice because R is also constrained by the self-gravity 
of the star. That is, we must minimize CAot = Uf + U (see eq. 
[3| with respect to R, not U alone. Generally, Uf cannot be a 
monotonic function of R, otherwise no stable unmagnetized 
hydrostatic configuration would exist (which is a necessary 
a priori assumption) . To illustrate this procedure, let us as- 
sume that U <^Uf, and thus the correction to R associated 
with the magnetic pressure is small (which will be shown 
below). Then, we may approximate Uf and U by 

Uf ~Uf,o + \^^5R^ and (7~?7o-[/o^, (36) 



Ro 



where 5R is the perturbation in the stellar radius as a result 
of the presence of the magnetic field. Thus, variations with 
respect to 5R give 

Uo/Rq Uo 



SR'. 



-Ro- 



(37) 



d^Uf/dR^ Uffi 

Therefore, subdominant fields will correspond to a small per- 
turbation to _Ro while equipartition fields will make correc- 
tions of order unity. However, it is clear that a stable config- 
uration with finite R and non- vanishing magnetic field does 
exist, and may be described by the analysis presented here. 



10 - 




Figure 2. The energy per unit helicity as a function of the force- 
free constant, cx, both appropriately scaled by the stellar radius. 
The solid, short-dashed and long-dashed lines represent the first 
three multipoles, Z = 1, 2, 3, respectively. The dipole configuration 
with aR = IT, labelled (a), corresponds to the global minimum 
energy for a given helicity; the associated field geometry is shown 
in Figure [T^. There is a spectrum of higher energy local minima 
for each multipole, e.g., (b) and (c) for the dipole, which corre- 
spond to Figures[Tj3 and[TJ;. The point labeled (d) is the minimum 
energy quadrupole configuration, and corresponds to Figure [TJl. 
For reference, the diagonal line corresponds to U = aH . 



3.4 Non-Equilibrium Configurations 

While we have found a globally optimized magnetic field 
geometry for a given magnetic helicity, thus far we have 
presumed that the star will naturally seek this minimum en- 
ergy state via large scale motions and reconnection events. 
However, there is no guarantee that the star will not be- 
come trapped in higher energy, local equilibria, in which the 
necessary large scale rearrangements are inhibited. For this 
reason we consider a simple evolution, comprised of a single 
multipole and varying a, which may occur, for instance, in 
some forms of current decay. 

For values of a which are not zeros of ji{aR) the ac- 
tion does not vanish. Nonetheless, it is straightforward to 
compute U and use S{aR) to relate this to H{a): 



U 



H 



U -S 



1(1 + 1) 



2 2 
Q a; 



2 

aai 



21 ji-i{aR) jtijaR) 
aR ji{aR) jf{aR) 



1 - 



2l + lji_i{aR) , jt^aR) 



+ 



aR ji{aR) jf{aR) 



,(38) 



(39) 



Since equilibria minimize the energy per unit helicity, we 
show in Figure[2]the quantity RU /H as a function of aR for 
a number of multipoles. For each multiple there is a global 
minimum, with U/H being smallest for the dipole. However, 
there are also a number of local minima corresponding to 



higher energy states. The explicit geometry for some of these 
are shown in Figure 1, and their implications are discussed 
in ED and ^ 



3.5 Solutions With no Toroidal Surface Stress 

The solutions discussed so far are true minima of the total 
energy subject to a fixed helicity. Yet, they are not true 
equilibria of the system since they have unbalanced stresses 
at the stellar surface. This is a paradoxical situation — how 
can a minimum energy configuration not be in equilibrium? 
The explanation is that the motions that the surface stress 
will induce motions that cause an unwinding of the toroidal 
component of the magnetic field and thus a reduction in the 
magnetic helicity of the system. In effect, the unwinding of 
the field causes helicity to flow out of the star to infinity. 
Thus, while the solutions discussed so far are indeed true 
minimum energy states so long as helicity is conserved, they 
are not relevant for a fully fluid star because such a star does 
not preserve helicity. 

One way to ensure helicity conservation is to make sure 
that there is no surface tangential stress, i.e., the toroidal 
magnetic field in the stellar interior vanishes as r ~* R. 
By equation (|26|l . this is easily arranged by having a field 
configuration with a single multipole I and choosing a such 
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/ 

/ 

/ 

- / 
/ 

^ 

5 10 

aR 

Figure 4. The same as Figure (2] except the zcro-torquc con- 
figurations shown in Figure |3] are labeled (e), (f), (g) (the first 
three dipole zero-torque states) and (h) (the first quadrupole 
zero-torque state). For reference, the diagonal line corresponds 
to [/ = olH. Note that these solutions are not at minima in U/ H , 
and thus are likely unstable. 



that 

ji{aR)=Q, ai^Q. (40) 

The situation is very similar to that of the minimum energy 
solutions discussed earlier, except that we now make use of a 
root of a different spherical Bessel function (cf. egs. 1341 and 
I40p . These are simply the standard spheromak solutions, 
though we shall refer to them as zero-torque solutions to 
contrast them with the other solutions discussed throughout 
this paper. 

As before, we have the option of choosing any value of 
I and any zero of the corresponding Bessel function. Thus, 
we have an infinite number of possible solutions, of which 
that corresponding to Z = 1 and the first zero of ji [aR) rep- 
resents the overall lowest energy state. All these solutions 
have vanishing field in the exterior vacuum, as required by 
the matching condition (|27p . and they have no unbalanced 
tangential stress at the surface of the star. However, these 
configurations still require the presence of a rigid bounding 
crust. This is because they do have an unbalanced radial 
stress resulting from the discontinuity in the tangential field 
at the surface (it is non-zero inside and vanishes outside). 
Thus these configurations are also inappropriate for an en- 
tirely fluid star. 

Figure [3] shows meridional slices of the magnetic field 
corresponding to the first three zero-torque solutions with 
dipole geometry (/ = 1) and the first quadrupole solution 
{I = 2). These solutions are very similar to the minimum 
energy solutions shown in Figure 1. The main difference is 
that the magnetic field is contained entirely inside the star 
and the field outside the star vanishes. Figure|4]is analogous 



Table 1. Approximate values of aR associated with various mul- 
tipole solutions for the minimum energy (ME) and zero-torque 
(ZT) solutions, described by equations II41I I and II43I I. Values asso- 
ciated with images shown are labeled by the corresponding figure 
number in the Notes column. 



Model Type n I aR Notes 






1 


7r 


Global Energy Minin 


1 


1 


2-K 


Fig.[Tj3 


2 


1 


37r 


Fig.IIJ: 





2 


4.4934 


Fig.[lji 


1 


2 


7.7253 




2 


2 


10.904 







1 


4.4934 


Fig.[3]3 


1 


1 


7.7253 


Fig.E 


2 


1 


10.904 


Fig. [3k 





2 


5.7635 


Fig.[3h 


1 


2 


9.0950 




2 


2 


12.323 





, Fig.[TJi 



to Figure 2, except that it highlights the solutions corre- 
sponding to the four panels in Figure O We note two inter- 
esting facts. First, the new solutions still satisfy U = aH; 
hence, for a given H , the value of a directly gives the energy 
of the configuration. The lowest energy solution is the one 
with the smallest value of a, i.e., solution (e) in Figures [3] 
and[4l Second, we see that the solutions do not live at min- 
ima of the energy with respect to a but rather at inflection 
points. This is because we determined the value of a from 
a boundary condition aX r = R rather than by minimizing 
the energy. 

A system that reaches one of these zero-torque solutions 
will conserve helicity since there is no outward flow of he- 
licity through the surface. However, it is not clear how long 
the system will survive in this state. Since there are neigh- 
boring solutions (with other values of a) which have smaller 
values of the energy, the system is likely to evolve away from 
the zero-torque solution. Once it does this, the stress at the 
surface will no longer vanish and helicity will cease to be 
conserved. Whether or not this sequence of events can hap- 
pen depends on the nature of the dynamical constraints that 
the fluid must satisfy, a topic that is beyond the scope of this 
paper. 



3.6 Explicit Expressions for the Magnetic Field 

In this subsection we provide for convenience explicit ex- 
pressions for the magnetic field geometries of a few low-order 
multipole solutions, including those presented in Figures[TJi,- 
d and[3j3-f. 

A natural normalization of the expressions obtained in 
sections 13. II is in terms of the standard multipole field com- 
ponents of the surface magnetic field, Bi. In this case the 
exterior and interior field of a single multipole component is 
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Figure 3. Meridional slices of illustrative magnetic field geometries with vanishing surface stresses for (e) the lowest energy such 
configuration, corresponding to the first dipole state, (f ) & (g) the second and third dipole solutions, and (h) the lowest energy quadrupole 
configuration. The poloidal and toroidal field structure are shown by the field lines and fiUed-color contours, respectively. In all plots the 
field line density is proportional to the poloidal fiux densities and the color scheme is proportional to the toroidal fiux densities. Both 
are normalized such that all frames have identical total helicities. Note that there is no magnetic field outside the star. 



2;iven by 

B,.>i{ — Bi 



R 



1+2 



Y, - 



l + l \r 



R 



1+2 



rji{aR) 

aRrji-i{ar) ~ lRji{ar) 



l{l + l)rji{aR) 

aR ji(ar) _ 
" l{l + l)ji{aR)' 



(41) 



The field may then be explicitly constructed given a value of 
aR, and the explicit forms of the vector spherical harmonics 
and spherical Bessel functions given in Tables lEll and IE2I 
respectively. The values of aR corresponding to the first 
few minimum energy (ME) multipole solutions are given in 
Table [U 

However, the normalization in equation pif) is diver- 
gent when the radial component of the surface field vanishes, 
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and thus for the zero-torque solutions described in section 
13.51 and shown in Figures [3] & 131 For these a better choice 
of normahzation is in terms of the contribution of a given 
multipole to the total magnetic helicity, which simplifies to 



10" 



aR 



(42) 



and also diverges for the zero-torque configurations 
[ji{aR) = 0] when the strength of the radial component 
of the surface field (Bi) is held fixed. Explicitly, in terms of 
the contribution to the helicity, the solutions with vanishing 
surface stresses are given by 



Br 



= 0, 



E 



'-Hf 



Rji jar) 
rji-i{aR) 



+ 



aRrji-i{ar) — lRji{ar) 



i{aR) 
aR ji (ar) _ 
1(1 + 1) ji-i{aR)' 



(43) 

Note that since we used the vanishing surface stress limit for 
the helicity, the above expression is only valid for the config- 
urations discussed in section 1X51 A more general expression 
can be normalized in terms of the helicity [eq. (|39p ] , however 
this is correspondingly more complicated. The values of aR 
associated with the first few multipoles of the zero-torque 
(ZT) configurations are given in Table [T] 



3.7 Star With a Rigid Crust 

We now consider a star with a rigid crust that is able to 
support unbalanced magnetic surface stresses. The partic- 
ular astrophysical system we have in mind is a magnetized 
neutron star. A rigid crust provides a firm anchor for any 
field lines that penetrate it. As a result, kinks in the mag- 
netic field at the surface can survive indefinitely and there is 
no tendency for helicity to flow out to infinity. An additional 
effect of the crust is that any radial field component that is 
present at the stellar surface when the crust first forms will 
survive unchanged for a long time. This permanent surface 
field serves as a boundary condition at r = i?. 

We are then led to pose the following new problem: 
Given a fixed structure of the radial magnetic field at the 
stellar surface (i.e., values of the coefficients ai) and a fixed 
total helicity of the system (the value of H), what are the 
equilibrium configurations of the magnetic field? 

The magnetic field distribution exterior to the star is 
uniquely determined by the coefficients ai. The solution in 
the interior of the star must be force-free and must there- 
fore take the form (|26|) . with the coefficients di given by the 
matching condition (|27|l . Therefore, the only free parameter 
in the solution is a, and its value must be chosen such that 
the system has the required helicity H. The problem is thus 
quite similar to that posed in i|3.5l The only difi'erence is 
that, instead of having a vanishing field at the stellar sur- 
face, we now have a prescribed non- vanishing distribution of 
field strength at the surface. 

It is useful to consider a specific numerical example. In 
the following, we have chosen the ai to scale as + 




Figure 5. The helicity as a function of aR in units of the 
B'iR*' /Sit where B* is the normalization of the surface dipole 
field. The thin sohd Unes are the helicities of the individual 
modes, normalized by the Kolmogorov scaling described in the 
text. The thick solid line is the total helicity associated with the 
Kolmogorov spectrum. Also shown is a constant helicity of 100 by 
the long-dashed line, which crosses the helicity curve at a num- 
ber of points. Note that the spikes in the helicity are due to the 
fact that we have fixed the magnitude of the exterior field, which 
implies that the interior field strength diverges at the zero-torque 
configurations. 



and thus Sf oc , to crudely model the effect of a Kol- 
mogorov turbulent cascade. Given a set of ai and a selected 
value of a, equation (|39p allows us to calculate the helicity 
of each multipole and equation (|38|l gives the correspond- 
ing energy. Moreover, since the vector spherical harmonics 
are mutually orthogonal, the total helicity and energy of the 
system are obtained by simply summing the contributions 
from the individual multipoles. 

Figure O shows the variation of total helicity (in appro- 
priate units) as a function of aR for the numerical example 
described above. The many spikes in the curve correspond 
to the various zeros of ji{aR); as equation (|39|) shows, the 
helicity diverges at these zeros. If the system has a certain 
conserved helicity H, which we have chosen to be 100 in 
the particular example shown in Figure O there will in gen- 
eral be multiple solutions for a which satisfy this constraint. 
Each of these is a valid solution to the problem. The solu- 
tions come in pairs, each bracketing a particular zero of a 
particular spherical Bessel function. 

Figure [51 shows the energies of these solutions as given 
by equation p8|l . We see that the solutions are close to, and 
bracket, the zero-torque solutions shown in Figure (H There 
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Figure 7. Meridional slices of a particular realization of the Kolmogorov spectrum described in the text for the states marked in Figure 
[B] each corresponding to a heUcity of 100. In particular, (i), (j) and (k) show the first, third and seventh lowest energy states with 
a helicity of 100, while (1) shows the more complicated structures that can appear in higher energy states. The poloidal and toroidal 
field structure are shown by the field lines and fiUed-color contours, respectively. In all plots the field line density is proportional to the 
poloidal flux densities and the color scheme is proportional to the toroidal flux densities. Both are normalized such that all frames have 
identical total helicities. Note the similarities with Figure [3] Note also that all the panels have identical exterior field distributions. 



is thus a close connection between the rigid-crust solutions 
and the zero-torque solutions. Finally, Figure[7]shows merid- 
ional slices of the magnetic field corresponding to the four 
solutions identified as (i), (j), (k) and (1) in Figure [5] The 
family resemblance of these solutions to those shown in Fig- 
ure [S] and indeed also Figure [TJ is obvious. 

The strong similarity between the rigid-crust and zero- 



torque solutions is easy to understand. We have selected a 
somewhat large value of H for the former, which means that 
the solutions must have a significantly larger field strength in 
the interior of the star compared to the surface, as is evident 
in Figure [T] In this situation, holding the surface field fixed 
at a small value (rigid-crust solutions) is nearly the same as 
setting the surface field to zero (zero-torque solutions). 
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Figure 6. U/H as a function of aR for the fixed Kolmogorov 
surface field spectrum described in the text. The points are states 
with H = 100 as shown in Figure [5] and those that are expUcitly 
shown in Figure[7]are labeled (i), (j), (k) and (1). Note that these 
are not necessarily energy minima, which is a direct result of fixing 
the multipolar structure at the stellar crust. 



creases, but the structure (described by a) remains un- 
changed. 

However, if the poloidal currents decay more rapidly 
than the toroidal currents, e.g., due to enhanced dissipation 
of the surface currents or because dissipation is non-linear, 
Q will be time dependent. This is clear from equation (|26|l . 
m which B"^ oc a. Thus if poloidal currents (and hence 
toroidal magnetic fields) decay more rapidly than toroidal 
currents (and hence poloidal magnetic fields) , d < and the 
magnetic field configuration will move to the left in Figures 
(2] 4 and 6. When the system evolves out of a local equi- 
librium state a substantial amount of magnetic energy is 
necessarily liberated. The natural timescale over which this 
energy is released is the Alfen crossing time, though we do 
not address the details of the instability here. 

Finally, note that in the presence of a conducting mag- 
netosphere with an infinitely conductive surface a fiuid star 
would appear as a uniform medium. This is despite the fact 
that the magnetosphere has a much lower density than the 
stellar interior, and thus is always described by force-free 
dynamics. This disparity results in very different dynam- 
ical times for the two regions. Nevertheless, since currents 
can freely flow between the stellar interior and the magneto- 
sphere, the minimum energy equilibrium configuration will 
be described by that for a uniform medium, which has no 
stable equilibrium magnetic field configuration for any value 
of magnetic helicity. That is, in this case a = U /H can take 
any value, including zero. Thus such a field is unconditional 
unstable. Whether or not this occurs in the context of a 
stellar field depends upon the conductivity of the surface, 
the exterior, and the ability of the crust to fix the magnetic 
field. 



3.8 Long Term Field Evolution 

In the presence of a non-zero resistivity, the magnetic field 
will necessarily decay. Since, in practice, the resistive de- 
cay timescales are typically much longer than the Alfven 
crossing time, we may approximate the state of the mag- 
netic field as a series of non-equilibrium configurations, for 
which the general solution is given by equations (I26|) and 
H2ip . Furthermore, we will assume that the field begins in a 
local minimum energy configuration, where the large-scale 
structure prevents further fast reconnection. 

As the field evolves due to resistive decay, the energy 
per unit helicity, U /H, will change as well. As a consequence, 
the value of a will evolve. This will generally result in the 
evolution of the magnetic field structure, provided the resis- 
tivity is not scale-invariant, i.e., small scale currents damp 
more rapidly than large scale currents. To see this consider 
Ohmic decay with a constant resistivity, p. In this case, the 
local field strength evolves according to 

2 

B = -VxE = -pVxJ = -AvxVxB = -^B. 

4-K Air 

(44) 

Upon inserting equation (|26|l into this, we find 

ill ~ — -. — ai and a = , (45) 
47r 

i.e., the magnetic field strength (proportional to ai) de- 



4 DISCUSSION AND APPLICATION TO 
NEUTRON STARS 

4.1 General Discussion 

Although we have described three different kinds of solutions 
— minimum- energy solutions ( i]3.2|l . zero-torque solutions 
(333} and rigid-crust solutions f i]3.7[) — all have certain 
features in common. 

For each class of solutions, the lowest energy state 
("ground state"), subject to a fixed total helicity, is dom- 
inated by a dipolar field configuration. Nevertheless, within 
the context of a magnetic relaxation model, the presence of 
higher order multipoles is not unexpected. After the cessa- 
tion of a dynamo in the early history of a star, the result- 
ing magnetic field will be chaotic and highly tangled over 
a wide range of scales. The field will also be violently un- 
stable, decaying rapidly towards a local minimum energy 
configuration. However, as described in tj3.8l this will only 
proceed as long as the rearrangement is driven by small scale 
reconnection. Large scale rearrangements, such as convert- 
ing low-order multipoles into dipoles, will occur over the 
much longer resistive timescale. As a consequence, we gen- 
erally expect the resulting field to be dominated by a dipole 
component with significant quadrupolar and octupolar com- 
ponents. Th is is consistent with the ou tcome of numerical 
simulations l|Braithwaite &: Spruitll2004l ). 

A novel feature of our work is the identification of 
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higher-energy quasi-equilibrium states. At late times, small- 
scale current decay will drive the system to evolve towards 
smaller a, necessarily resulting in periods of violent instabil- 
ity when the system hops from one quasi-equilibrium state 
to the next. These will result in large scale reorganizations of 
the internal and external stellar magnetic fields over Alfvenic 
timescales. As a consequence, these events could power ac- 
tive periods in which a substantial fraction of the stellar 
magnetic energy is released. 

Since the quasi-equilibria of the dipole state are located 
at zeros of low-order spherical Bessel functions, the magnetic 
energy release in these events should be roughly quantized 
(though the system may jump over many local equilibria in 
a single event), releasing an energy of order kivH/R for some 
integer k. The observed energy release will depend strongly 
upon the emission mechanism and will likely not appear 
strictly quantized. Nevertheless, this model implies a lower 
limit upon the energy that can be released during an out- 
burst, simply given by the energy of the global minimum 
energy configuration. 

The minimum-energy solutions described in i]3.2l and 
Figures 1, 2, have a toroidal kink in the magnetic field at 
the stellar surface, in apparent contrast to numerical simu - 
lations fcf. lYoshida et al.ll2006l :[ Braithwaite fc Spruitll2004 ). 
The toroidal field discontinuity will lead to a non-radial sur- 
face stress, which cannot be balanced by a fiuid star. Thus, 
these solutions are probably not of practical interest. 

The zero-torque solutions described in i)3.5l and Figures 
[31 131 are more promising. They are perfectly force- free solu- 
tions with no unbalanced stresses. However, they have some 
deficiencies. First, they have no magnetic field outside the 
star. Therefore, they are not very useful for understanding 
the observed exterior field of systems such as magnetic A 
stars and white dwarfs. Second, these solutions are not min- 
imum energy states of the force-free problem. Therefore, a 
star in one of these states could, in principle, lower its energy 
by evolving to a nearby configuration with a lower value of 
a. The neighboring lower-energy state will not be stress-free 
at the surface and the system may thus develop a dynam- 
ical instability. Whether or not this will happen in a real 
system depends on dynamical constraints which we cannot 
address. Modulo this caveat, the zero-torque solutions are 
of possible interest for understanding the numerical work of 
iBraithwaite fc Spruit! (|2004l ). 

All the solutions we have described in this paper have 
a constant a in the stellar interior. This is a natural con- 
sequence of the variational approach we have taken, which 
is ba sed on the principle of conservation of total helicity 
ifravlor .igT^ l. but it is not necessarily appropriate for a 
magnetized star. We could, in principle, consider solutions 
in which a part of the stellar interior has a non-zero value of 
a and contains all the toroidal field and helicity, while the 
rest of the star has a = and purely poloidal structure. Such 
solutions are more difficut to construct, but may be able to 
produce stable configurations with no stress at the surface. 
An even more extreme scenario is to consider a separately 
conserved microscopic helicity for each magnetic fiux tube 
(see Appendix A) . These topics are beyond the scope of this 
paper. 



4.2 Neutron Stars 

During their birth, neutron stars are expected to have con- 
vectively driven dynamos and to develop extremely strong 
magnetic fields, as high as 10^^~^^ G. However, shortly af- 
terward, and prior to the development of a solid crust, the 
dynamo action ceases and the nascent neutron star magnetic 
field evolves for roughly 100 s. During this time magnetic re- 
laxation will play a central role in determining the strength 
and structure of the magnetic field. Once the crust forms, the 
magnetic field at the surface of the star is frozen in and the 
magnetic helicity of the system is fixed. Further evolution of 
the system is likely governed by total helicity conservation, 
coupled with the requirement that the radial field at the 
stellar surface is fixed. These are precisely the conditions we 
have assumed in deriving the rigid-crust solutions discussed 
in i]3.7l and Figures 5-7. 

Generally, we expect a dominant dipole field structure, 
with non-negligible quadrupole and octupole components. 
Since the magnetic helicity is a scalar quantity, we naturally 
expect a non-aligned magnetic field with respect to the rota- 
tion axis. For a highly tangled initial field (with considerable 
power on small-scales), we would expect the orientation of 
the birth fields of neutron stars to be uniformly distributed 
relative to the spin axis. This provides a natural way in 
which to produce the off-axis fields required for understand- 
ing radio and X-ray pulsars. 

For a system with non-zero total helicity, since any mag- 
netic field at the stellar surface will generally have an asso- 
ciated toroidal stress, we speculate that the fiuid neutron 
star before the formation of the crust will evolve through 
states in which the surface field is generally weak relative 
to the field in the interior. This condition will be frozen 
in once the crust forms. Thus, the frozen-in helicity will be 
large compared to the characteristic helicity we might expect 
based on the magnitude of Bs, the surface field strength: 
H 3> (-Bg /87r)_R*. Similarly, the magnetic energy will also 
be large: U > {B^/8tv)R^. These are the conditions we as- 
sumed for the examples shown in Figures 5-7. 

Once the crust forms, the system will settle down into 
the nearest available solution, which would presumably cor- 
respond to a relatively large value of a. Subsequently, we 
imagine that the magnetic field will evolve via sudden jumps 
from one solution to the next, with a steady reduction in the 
magnitude of a. Each jump will cause a decrease in the total 
magnetic energy and a corresponding increase in the thermal 
energy of the star. The leaking of this thermal energy out 
of the star will power thermal X-ray emission from the sur- 
face. This may explain the quasi-steady X-ray emission seen 
in anomalous X-ray pulsars (AXPs) and other magnetars. 

It is interesting to note that the source of this thermal 
energy is by its very nature non-steady. Bursts of energy are 
released each time the interior magnetic field switches from 
one state to the next. Although the fiux escaping from the 
surface will be smoothed on the thermal diffusion time, some 
variability is expected and is indeed seen in magnetars. 

Furthermore, as we see from Figures 5 and 6, neighbor- 
ing solutions are bunched more closely together for larger 
values of a and are farther apart at lower a. Since the mag- 
netic energy is to a good approximation proportional to a 
(see Fig. 6) , we expect frequent small bursts of energy early 
in the life of the neutron star, when a is large, and rarer 
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but more powerful bursts later on. Thus, there should be 
less variability in younger systems and larger fluctuations in 
older systems, with the most dramatic variations occurring 
in the oldest magnetars. 

The above discussion was concerned with magnetic field 
evolution in the interior of the star with a fixed surface field. 
However, since the star has a much sronger field in its inte- 
rior compared to the surface, we expect the crust occasion- 
ally to crack under the pressure and to release some of the 
enclosed energy and helicity in a sudden cataclysmic event. 
These events dump a lot of energy near the surface of the 
star and provide a natural mechanism for powering bursts 
associated with soft-gamma repeaters (SGRs). With typical 
energies ranging from io^®~*® erg, SGR bursts require inte- 
rior magnetic fields on the order of 10^~^* G, which are con- 
siderably less than the equipartition values that are possi- 
ble a s a result of m agnet ar formation (jThompson fc DuncanI 
ll993l : ISpruij|2002D . For magnetars, the Alfven crossing time 
is on the order of a second, corresponding to maximum lumi- 
nosities on the order of 10^^ erg s~^ . If this energy is thermal- 
ized in the neutron star crust over the Alfven timescale, it 
would result in surface temperatures on the order of roughly 
an MeV, similar to what is seen in SGR fiares. However, the 
precise details of the emission will d epend upon the details 
of th e emission mechanism (see, e.g.. [Thompson et al.|[20o3 : 
iLvutikov 2006 ). 

As in the case of internal field rearrangements, we ex- 
pect the most violent events to occur somewhat late in the 
life of a magnetar, when the interior field is dominated by 
lower multipoles. Thus, the strongest outbursts should be 
seen in relatively old systems. All of these arguments sug- 
gest that SGRs as a class might be older than AXPs. 

In this picture magnetars are born with magnetic fields 
in configurations that are local equilibria at strengths orders 
of magnitude larger than the global minimum. If we assume 
that the minimum energy of observed SGR bursts is indica- 
tive of the energy scale of the power source, this implies that 
the field strength in this minimum energy state is roughly 
10^ G. Hence, the SGR bursts would be the result of the de- 
cay of the internal magnetic field, resulting finally in a stable 
field strength commensurate with millisecond pulsars. Note 
that this is a different mec hanism than that described in 
iBraithwaite fc Spruit! (|2006l ). 
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APPENDIX A: RELATIONSHIP TO OTHER 
VARIATIONAL PRINCIPLES 

We present a particular variational principle, with the some- 
what surpris ing result (though similar to that found in 
IWoltieij ll962l') that the minimum energy states of the mag- 
netic field correspond to a completely force-free state in 



the stellar interior with uniform force-free constant, a. This 
is surprising since equilibrium configurations that balance 
magnetic stresses against fiuid pressure are conspicuously 
absent from our approach. Thus it is worthwhile to discuss 
the particular conditions that our solution is valid under and 
why this differs from those derived under the assumption of 
ideal MHD. 

Our v ariational p rinciple differs explicitly from that de- 
scribed in lFieldlll986l (an interesting paper brought to our 
attention by the anonymous referee) . The latter reproduces 
the standard MHD relation, balancing pressure against mag- 
ne tic stresses, lea ding the author to argue that the result 
of IWoltied l| 196211 is only vahd if the pressure is uniform 
throughout the fluid. This variational principle is derived 
in a fashion nearly identical to that presented in section [2j 
minimizing the total energy (fiuid and magnetic) subject to 
a condition upon the helicity. However there are two impor- 
tant distinctions. 

The fir st concerns t he terms in the fiuid component of 
the energy. iFieldl (|l986') ignores the contributions due to 
self-gravity. As a consequence pressure gradients can only 
be balanced by the Lorentz force. This generally will result 
in an unconditionally unstable system for isolated magne- 
tized fiuids. In contrast, a self-gravitating magnetized fiuid 
can balance pressure gradients against the gravitational ac- 
celeration, and thus the Lorentz force may in principle van- 
ish. Therefore, when we include self-gravity, the existence of 
pressure gradients is not synonymous with the magnetic field 
violating force-free (see, e.g., [Chandrasckhar & Prcndcrgas^ 
119561 ) 

The second concerns the way in which helicity conserva- 
tion is imposed. In ideal MHD it is possible to define a con- 
served microscopic magnetic helicity associated with each 
magnetic fiux tube. Therefore, we may in principle seek to 
impose conservation of not only the macroscopic, volume in- 
tegrated helicity, but the helicity associated with each field 
loop separately, i.e. take a not to be a constant, but to be a 
fu nction of pos ition. This is precisely the constraint imposed 
in lFieldl (|l986l l. 

However, it is far from clear that in the presence of 
large-scale dissipation that this is the co rrect prescription. 
Indeed, the conjecture by I Taylor! (|l974 ). apparently sup- 
ported by spheromak experiments, is that this be replaced 
with a volume-integrated con straint, of the form we have 
employed. The reason given bv lTavloil (Il974l ) for this is that 
the microscopic helicity associated with each fiux tube is not 
conserved during reconnection events. Indeed, during recon- 
nection, fiux tubes with different helicities may be combined, 
changing the total helicity of the new loop (the precise man- 
ner in which this occurs is not clearly defined since the evo- 
lutionary equation for the fiux-tube helicity depends upon 
the way in which helicity proceeds). However, reconnection 
will not change the macroscopic topology of the magnetic 
field, and therefore the total helicity will be preserved. 

An alternative way to justify this is to note that if the 
magnetic field rapidly reconnects, it is impossible to define 
unique magnetic field loops. That is, after sufficient time 
each magnetic field loop will have been connected with ev- 
ery other remaining loop at some point in the reconnection 
history. As a consequence, the force-free constant (i.e., a) 
associated with each loop will necessarily be the same, since 
this is constant along magnetic field lines. This will remain 
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to be the case unless the hehcity along some subset of field 
lines can preferentially decay (which is the case in the ab- 
sence of a rigid crust). 

These differences explain the absence of solutions which 
balance magnetic stresses and pressure gradients in our for- 
malism. Such solutions require additional constraints, such 
as flux freezing in ideal MHD, which inhibit the movement 
of magnetic field lines through the fluid. By considering 
the high-dissipation/rapid-reconnection limit via our vol- 
ume constraint upon the total magnetic helicity we have 
necessarily ignored such constraints. That is, in the absence 
of such additional constraints there are no stable configura- 
tions which balance magnetic stresses against pressure gra- 
dients, since the fluid and magnetic field will move through 
each other until such gradients are erased. (This is identi- 
cal to equilibrium configuration of a two-component ideal 
gas, in which both species separately satisfy the equilib- 
rium condition.) In the rapidly reconnecting limit, the mag- 
netic field can globally reconfigure on short timescales, and 
thus we are justified in treating the magnetic field and fluid 
as non-interacting in this fashion. Finally, because we are 
restricting our attention to gravitationally bound objects, 
non-vanishing pressure gradients may be balanced against 
self-gravity, and therefore do not require violations of the 
force-free condition to maintain equilibrium. 



APPENDIX B: FORCE-FREE FIELDS IN 
UNIFORM MEDIA 

When the medium is uniformly infinitely conductive (i.e., 
the currents required by the resulting magnetic field satisfy 
any physical restrictions) the quantity to be minimized is 



Stt 



(B ■ B - aA ■ B) dx , 



(Bl) 



where the constraint upon the total helicity is incl uded via 
the La grange multiplier a. As originally discussed in lWoltieil 
l) 19581 ) (and thus frequently referred to as the Woltjer state) 
variations with respect to A produce the familiar force-free 
equation VxB = qB for the magnetic field. That this pro- 
duces a force-free field is not unexpected; non-force-free per- 
turbations to a force-free solution necessarily require work 
to be done upon, and thus energy added to, the magnetic 
field. For the purpose of comparison with what follows we 
will derive this equation below. 
First note that 



5{B ■ B) = 25A ■ (VxB) - 2V-(B x SA) , 
5(A • B) = 25A ■ (VxA) - V-(A x SA) . 

Therefore, 

8tt5S ^ I 5A- (VxB-aVxA)dx 
Jv 



(B2) 



BxSA AxSA] -dS 

2 



= [ SA - (VxB-QB)d^a; 
Jv 

- (f SA-[. 

Jav L 



(B3) 



fix (B-|A)|d\, 



is generally gauge dependent. This can be cured if SA is 
constrained to vanish at dV . Since in what follows we will set 
dV to infinity and the energy in the field will necessarily be 
finite, a gauge will necessarily exist in which A vanishes on 
dV , and thus this constraint upon dV is justified. Therefore, 
the desired magnetic field must satisfy 

VxB = aB. (B4) 

Equation HB4|I implies that 

Vx(B-aA)=0 ^ aA^B-faVA, (B5) 

and thus a natural gauge choice for non-zero a is the "force- 
free" gauge defined by A = 0. Note that this is also a Lorentz 
gauge since V-B = 0. 

Within the force-free gauge it is straightforward to show 
that the helicity and minimum possible magnetic energy are 
proportional. That is, 

H = [ — A Bdx= / — Q~^B -Bd^a; = a"^C/, (B6) 
Jv Jv 

and thus small a correspond to small U for the same H. 



APPENDIX C: BOUNDARY CONDITIONS 

On the boundaries between type-I and type-II regions, the 
fields must satisfy the usual boundary conditions: 

47r 

n-(Bi-B„)=0 and n x (Bi - B„) = — K , (CI) 

c 

where n is the surface normal (directed from I to II) and K 
is the induced surface current density. The proper surface 
current can be found by explicitly constructing solutions and 
minimizing the magnetic energy. Note that this ignores any 
energy associated with the surface currents themselves. 

However, this should also emerge naturally from the 
equations (|6]) and ((7]). The first is obtained from V-B = 0, 
which is true by definition and may be verified by taking the 
divergence of equation (|6}. The second is obtained via in- 
tegrating equation ([6]) around an infinitesimal loop crossing 
the surface. Define a local Cartesian coordinate system at 
the position of interest at the I-II boundary in which z || n 
at the boundary, at which z = Q. Then, upon integrating 
around a rectangular loop in the x-z plane of length I in the 
x-direction and vanishing length in the z-direction, we find 

(Bi -Bii) -^x = / B ■ d£ = / (VxB)-dS 

Jc J A 

= / (aB- VxVx/A) -dS 

J A 

= j [aA + A x V/ - /VxA) ■ dl 

[(1 - f)aA - /aVA + \xVf]-dl. 

(C2) 

The third term in the integrand can be integrating after 
noting V/ = i5(z) by definition, and thus 



5{z)\ X z ■ d£ = (A X z) ■ z = 0. 



(C3) 



where n is the normal to the surface defined by dV and 
d\ is defined by dS = nd^x. Note that the surface integral 



Despite the explicit presence of the vector potential, the first 
two terms in the integrand are indeed gauge independent. 
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Explicitly, note that since under the gauge transformation 
A ^ A + VA' also implies that A ^ A — A', and hence 

(1 - f)aA - /aVA ^ (1 - f)aA - faVA + VA' . (C4) 

Finally, noting that J^'V A' ■ d£ — completes the proof. 
Therefore, 



(Bi - Bii) ■ £x = qAi ■ est + a£x ■ VA 



(C5) 



With an analogous expression for a loop in the y-z plane, 
the surface current is given by 

47r . 



-if = an X (Ai + VA) 



(C6) 



If we choose the force-free gauge within regions of type-I, 
continuity of the vector potential gives 

47r . 



-K = nx (Bi +aVA) , 



fi X (Bii - aVA) =0, 



(C7) 



(C8) 



on the boundary. That is, as expected Bn is a potential field, 
though otherwise the surface currents are unconstrained. 



APPENDIX D: THE ACTION IN 
NON-UNIFORM MEDIA 

In non-uniform medium, it is possible to explicitly reduce 
the computation of the action to integrals over the bound- 
aries between conducting (type-I) and insulating (type-II) 
regions. To do this, first note that the magnetic energy is 
given by 

U = [ ^B-Bd^2:+ / ^i-B- Vi^d^a; 

1 (Dl) 
— ^B-dS, 

/e{i-ii) 

where V-B = was used. Note that because the n ■ B is 
continuous across the boundary, this may be evaluated using 
only type-I quantities. 

If Al is defined by A = Al + VAl where Al is in 
the Lorentz gauge (i.e., V-Al = 0), the magnetic helicity is 
given by 



Stt 



= / — B-Bd^i 

1 1 8-K 



Stt 



a 



I — B-Bd^2:- / 
J I 8^ Jg 



e{i-ii) 



-^A.d5 



—LpV-AA\ 
11 Svr 

-1 / _LB.Bd\-Q-^ / —ifi'B- AS 



Stt 



- / — v^V'ALd'x 



e(i-n) 



S(I-II) 



-!-V5VAl • dS / -^B-VA 



1 



f-1 

Jii Stt 



^a-^U-i -^(AlVv3"(/=VAl) -dS, 
ia(i-ii) fSTT 

(D2) 

where V-B = V-Al = and the continuity of A and B ■ n 
at the boundary were used (i.e., on the boundary n • An = 



n-Ai = a ^n-Bi = Q ^fi - Bu). Therefore, the action may 
be written as 

S = U~aH=i — (AlVv? - .^VAl) • dS , (D3) 
/e(i_ii) Stt 

where tp and Al are implicit functions of a through equation 
Note that when Al vanishes, i.e., the force- free gauge 
is a Lorentz gauge everywhere, U — aH as found for the 
Woltjer state. 



APPENDIX E: VECTOR SPHERICAL 
HARMONICS 

Vector spherical harmonics are used extensively, and thus 
some of their prop erties are summarized below (see, e.g., 
iBarrera et al.|[l985l , for more detail). While it is possible to 
expand the angular dependence of a vector in many ways, 
simplifications of the Helmholtz equation can be obtained 
with the following choices: 

Yin, = evYim , *im = rVYi™ , ^ira = rx VYi™ , (El) 

where the Yim are the standard scalar spherical harmonics 
(see Table lEll for explicit expressions for low I). These have 
the following useful properties: 



V-(F(r)*^ 



0, 



Vx {F(r)Yi^) = -^^i„,. 



Vx {F{r)^i„,) = ( -drrFir) ] 



(E2) 



Vx {F{r)^im) 



-F{r-) Y 



ia.rF(r) ) ^im . 



The are also orthogonal on the unit sphere: 



J dn*;„ ■ = l{l + l)5wS^rn' , 

= J dfl^i^ ■ = . 

Note that equation HE2|) implies that 



(E3) 



V x V x (F(r)*,^) = (-A,F(r))*^ 



(E4) 



where A; is the 3-dimensional Laplacian associated with the 
meridional harmonic I, i.e.. 



Note that the general solution to Aif + a^f = are 
the spherical Bessel functions of order I, ji{ar) and ni{ar). 
The first few of these are sum marized in Table |E2| and mor e 
information may be found in lAbramowitz fc StegunI l| 19721 ) . 
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Table El. Explicit expressions for the azimuthally symmetric 1 = 0,1,2 and 3 vector spherical harmonics. 




Table E2. Explicit expressions for the spherical Bessel functions of orders Z = 0, 1, 2 and 3 l|Abramowit2 fc Stegunlll972^ . 
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